An h-adaptive spacetime-discontinuous Galerkin method for linear elastodynamics

We present an h-adaptive version of the spacetime-discontinuous Galerkin (SDG) finite element method for linearized elastodynamics (Abedi et al., 2006). The adaptive version inherits key properties of the basic SDG formulation, including element-wise balance of linear and angular momentum, complexity that is linear in the number of elements and oscillationfree shock capturing. Unstructured spacetime grids allow simultaneous adaptation in space and time. A localized patch-by-patch solution process limits the cost of reanalysis when the error indicator calls for more refinement. Numerical examples demonstrate the method’s performance and shock-capturing capabilities. RÉSUMÉ. Nous présentons une version h-adaptative de la méthode de Galerkin discontinu espace temps (SDG) pour l’élastodynamique linéaire (Abedi et al., 2006). La version adaptative hérite des principales propriétés de la formulation de base de la méthode SDG, comme par exemple l’équilibre par élément de la quantité de mouvement et du moment cinétique, la complexité informatique linéairement proportionnelle au nombre d’éléments et la capture sans oscillation de chocs. Les maillages, qui sont non structurés en espace temps, permettent l’adaptation simultanée en espace et en temps. Une procédure de résolution localisée par patch permet de limiter le coût de la réanalyse après le raffinement local imposé par l’indicateur d’erreur. Des exemples numériques sont présentés et montrent l’efficacité de la méthode et plus particulièrement sa capacité à capter des chocs.